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Abstract 

We introduce a special stochastic perturbation of the flow of diffuse matter as a curve 
in the group of diffeomorphisms of flat n-dimensional torus such that the perturbed 
system yields a solution of Burgers equation in the tangent space at unit of the 
diffeomorphism group. The same perturbation of the flow of perfect incompressible 
fluid yields a solution of Reynolds type equation but under some special external 
force on the diffeomorphism group it transforms into a solution of Navier-Stokes 
equation without external force. 
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1 Introduction 

The paper is devoted to the Lagrangian approach to hydrodynamics in terms of geometry 
of groups of diffeomorphisms, suggested for perfect fluids by Arnold pQ and Ebin and 
Marsden |2] . In previous works by the author it was found that the adequate description 
of viscous fluids in this language requires involving stochastic processes such that their 
expectations are flows of viscous incompressible fluids (see, e.g., [31 H]). In this framework 
Newton's second law on the groups of diffeomorphisms, used in the case of perfect fluids, is 
replaced by its special stochastic analogue in terms of Nelson's backward mean derivatives. 
After transition to the tangent space at unit of diffeomorphisms group, there arises the 
Navier-Stokes equation via a natural modification of construction by Arnold, Ebin and 
Marsden that yields the Euler equation in the case of perfect incompressible fluid. In 
complete form this idea is realized in the case where the drift of above-mentioned process 
on the group of diffeomorphisms is right-invariant (see P,[7]). 

Here we consider another case. We introduce a special stochastic perturbation of a 
flow of diffuse matter, satisfying the above-mentioned stochastic Newton's law, and show 
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that the corresponding curve in the tangent space at unit satisfies Burgers equation. The 
same perturbation of a flow of perfect incompressible flow without external force satisfies 
the above-mentioned stochastic Newton's law as well, but yields a curve in the tangent 
space at unit that is a solution of a Reynolds type equation. Nevertheless, under the 
action of a certain special external force on the flow, this curve becomes a solution of 
Navier-Stokes equation without external force. We consider the fluid motion on the flat 
n-dimensional torus (we recall the definition below). 

The research was supported in part by RFBR Grants No. 07-01-00137 and 08-01- 
00155. 

This paper originated from a discussion with B. Rozovskii held long ago at the Uni- 
versity of Warwick. I am grateful to him for that initial prompt. 



2 Preliminaries 

Consider a stochastic process £(£) in M n , where t € [0,T], given on a certain probability 
space (n,J-,P) and such that £(£) is an L 1 -random variable for all t. The "present" 
("now") for £(£) is the least complete cr-subalgebra J\f£ of T that includes preimages of 
the Borel set of W 1 under the map £(£) : Q -> R n . We denote by £f the conditional 
expectation with respect to N%. The least complete a-subalgebra that includes preimages 
of the Borel set of lR n under all maps £(s) : f2 — ► M n for s < t (resp. s > t) is called the 
"past" (resp. "future") a-algebra and is denoted by V\ (resp. T\). 
Below we most often deal with the diffusion processes of the form 

+ / a{s^{s))ds + aw(t) (1) 
Jo 

in M n and flat torus T n as well as natural analogues of such processes on groups of 
diffeomorphisms (infinite-dimensional manifolds). In (pQ) w(t) is a Wiener process adapted 
to £(t), and a(t,x) is a vector field; a > is a real constant. 
Following Nelson (see, e.g., [8] - [10]) we give the next 

Definition 1 (i) The forward mean derivative D£(t) of the process £(£) at t is the L 1 - 
random variable of the form 



v ' 1 \ At ) 



(2) 



where the limit is supposed to exist in L 1 (f2, J 7 , P) and At — > +0 means that At — >• and 
At > 0. 

(ii) The backward mean derivative D*£(t) of £(£) at t is the ^-random variable 

^)^L^(^) (3) 

where (as well as in (i)) the limit is supposed to exist in L 1 (fi, J 7 , P) and At +0 means 
the same as in (i). 
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Notice that, generally speaking, D£(t) ^ D*£(t) (but, if almost surely (a.s.) has 
smooth sample trajectories, these derivatives evidently coincide). 

From the properties of conditional expectation it follows that D£(t) and D*£(t) can 
be represented as compositions of £(£) and Borel measurable vector fields 

^.)-jaj.*(ffi±MzM|« t ,-,), 

on M n (following [H] we call them the regressions): D£(t) = Y°(t,£(t)) and D*£(t) = 

Lemma 1 For a process of type (TJ]) ; we have D£(t) = a(t,£(t)) and so Y°(t,x) = a(t,x) 

For details of the proof, see, e.g., PE]. 

Let Z(t,x) be C 2 -smooth vector field on M. n . 



Definition 2 The L 1 -limits of the form 



DZitm) = Jon,* (*^™>) 

and 

D.Z( tl£W ) = ^ 

are called forward and backward, respectively, mean derivatives of Z along £(•) at time 
instant t. 



(5) 



(6) 



Certainly DZ(t,£(t)) and D^Z(t,^(t)) can be represented in terms of corresponding 
regressions, defined analogously to (J3J). If it does not yield a confusion, we denote these 
regressions by DZ and D*Z. 

Lemma 2 For the process <^ in M. n , the following formulae take place: 

DZ=^Z+(Y°-V)Z + ^-V 2 Z, (7) 

and 

D*Z =j t Z+ (*? ■ V)Z - y V 2 Z (8) 

where V = (gfr, gfsO> V 2 tae Laplacian, the dot denotes the inner product in W 1 , 
the vector fields Y°(t,x) and Y®(t,x) are introduced in Q). 

Recall that the n-dimensional torus T n can be considered as the quotient space of 
W 1 with respect to the integral lattice Z n . Introduce the Riemannian metric (•, •) on T n 
inherited form the Euclidean inner product in M. n . This metric is called flat and T n with 
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this metric is called the flat torus. Everywhere below we deal with the mean derivatives 
and fluid motion on flat torus. 

The main idea of description of viscous hydrodynamics in the language of mean deriva- 
tives is as follows. Consider the vector space Vect^ of all Sobolev i!P-vector fields 
(s > | + 1) on T n . 

Let a random flow £(t,m) with initial data £(0,m) = m G T n be given on a flat 
n-dimensional torus T n . Suppose that it is a general solution of a stochastic differential 
equation of the type 

£(£,m) = m+ / a(s,^(s,m))ds + aw(t). (9) 
Jo 

Let .D „,£(£, m) = v(t, m)), where v (t, m) is a C^-smooth in t and C 2 -smooth in m G T n 
vector field on T ra ; let a > be a real constant. Let £(t,m) satisfy the relation 

D*D^(t,m) = F(t,m), (10) 

where F(t, m) is a vector field on F n . Taking into account formula ([8]), we obtain 

D*D^(t, m) = (—v + (v, V)v - G — V 2 v) = —v + (v, V)v - yV 2 ^. 
Thus (flOl) means that v(t,m) satisfies the relation 

-v + (v,V)v-^-V 2 v = F(t,m), (11) 

that is the Burgers equation with viscosity ^- and external force F(t, m). We interpret fflUj) 
as a stochastic analogue of Newton's second law on the group of Sobolev diffeomorphisms 
V s {T n ) (see the next Section). 

The case of viscous incompressible fluids requires some additional constructions. In- 
troduce the L 2 -inner product in Vect^ by the formula 

(X,Y)= I (X(m),Y(m))fi(dm) (12) 

where (-, •) is the Riemannian metric on T n and fi is the form of the Riemannian volume. 

Denote by f3 the subspace of Vect^ consisting of all divergence-free vector fields. Then 
consider the orthogonal projection with respect to ([12]) : 

P : Vect {s) -> 0. (13) 

It follows from the Hodge decomposition that the kernel of P is the subspace consisting 
of all gradients. Thus, for any Y G Vect^ s \ we have 

P(Y) = Y - gradp, (14) 

where p is a certain iP +1 -function on T n , unique to within the additive constant for a 
given Y. 
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Let a random flow £(£, m) with initial data £(0, m) = m G T n be given on a flat 
n-dimensional torus T n . Let m) be the general solution of a stochastic differential 
equation of the type (j5J) and let D*£(t,m) = u(t,£(t,m)), where u(t,m) is a C 1 -smooth 
in t and C 2 -smooth in m G T n divergence-free vector field on T n , and let a > be a real 
constant. Suppose that m) satisfies the relation 

PD*D^(t,m) = F(t,m), (15) 

where F(t, m) is a divergence-free vector field on T n . Taking into account formulae © 
and f|T4|) . we obtain 

o 2 

PD*D4(t, m) = P(-u + (u, V)« - y V 2 «) 

= — it + (m, V)« — — V 2 u — gradp. (16) 

Thus (|T5l) means that w(t, m) is divergence-free and satisfies the relation 

d a 2 

—u+(u,V)u — —V 2 u — gr&dp = F(t,m), (17) 

that is the Navier-Stokes equation with viscosity ^- and external force F(t,m). 

We interpret f|T5|) as a stochastic analogue of Newton's second law on the group of 
Sobolev diffeomorphisms V s {T n ) of the torus, subjected to the mechanical constraint. 



3 Basic notions from the geometry of groups of dif- 
feomorphisms of flat torus 

The tangent bundle to T n is trivial: TF n = T n x M. n . Note that the flat metric generates 
in the second factor the inner product, same as in the copy of M. n from which the torus is 
obtained by factorization. 

Consider the set V s (T n ) of all diffeomorphisms of T n to itself belonging to the Sobolev 
space H s , where s > ~n + 1. Recall that for s > \n + 1, the maps belonging to H s class 
are C 1 -smooth. 

There is a structure of smooth (and separable) Hilbert manifold on V s (T n ) as well 
as the natural group structures with respect to composition. A detailed description of 
the structures and their interconnections can be found in [2]. Note that the tangent 
space T e T> s (T n ) at the unit e = id is Vect^ (see above). Recall that T e V s {T n ) contains 
its subspace (3 consisting of all divergence-free vector fields on T n belonging to H s (See 
Section [2]). 

The space TfT> s (T n ), where / G T> s (T n ), consists of the maps Y : T n — > TM such 
that nY(m) = f(m), where 7r : TT n — > T n is the natural projection. Obviously for any 
Y G T f V s {T n ) } there exists unique X G T e V s {T n ) such that Y = Xof. In any T f V s {T n ) 
we can define the L 2 -inner product by analogy with (fT2j) by the formula 

(X,Y) f = [ (X(m),Y(m)) f{m) »(dm). (18) 



5 



The family of these inner products form the weak Riemannian metric on V s (T n ) (it 
generates the topology of the functional space H° = L 2 , weaker than H s ). 

The right translation R f : V s {T n ) -> V s {T n ) } where R f (6) = 9 o / for 9, f G V s (T n ), 
is C^-smooth and thus one may consider right-invariant vector fields on T> s (T n ). Note 
that the tangent to the right translation takes the form TRfX — X o f for X G TT> s (T n ). 

The right-invariant vector field X on T) s (T n ) generated by a given vector X G 
T e V s (T n ) is C fc -smooth if and only if the vector field X on T n is # s+fc -smooth. This 
fact is a consequence of the so-called uj- lemma (see pj) and it is valid also for more com- 
plicated fields. For example, if a tensor (or any other) field on T n is C°°-smooth, the 
corresponding right-invariant field on V s (T n ) is C°°-smooth as well. 

Remark 1 The left translation L f : V s (T n ) -> V s (T n ), where L f (9) = 9 o f for 9,f G 
V s (T n ), is only continuous. Specify a vector x G M, n and denote by l x : T n — > T n the 
diffeomorphism l x (m) = m + x modulo factorization with respect to the integral lattice. 
Note that the left translation L\ x is C°° -smooth. 

Introduce the operators: 

B : TT n -> R n , 
the projection onto the second factor in T n x MJ 1 , and 

A(m) : R n — > T m T n , (19) 

the converse to 5 linear isomorphism of M n onto the tangent space to T n at m G T n . The 
map ^4 has the following property. For the natural orthonormal frame b in R n we have 
an orthonormal frame A m (b) in T m T n , the field of frames A(b) on T n consists of frames 
inherited from the constant frame b. Thus, for a fixed vector X G M n , the vector field 
A(X) on T n is constant (i.e., it is obtained from the constant vector field X on M. n and has 
constant coordinates with respect to A(b)) and, in particular, A(X) is C°°-smooth and 
divergence-free since such is the constant vector field X on M. n . So, A may be considered 
as a map A : W 1 -> /? C T e V s {T n ). 
Introduce 

Q 9(m) = A(<7(m)) o S, (20) 

where 51 G V s (T n ), m G T n . For a vector F G T f V s (T n ), we get Q ff y = A(p(m)) o 
B(y(m)) G T g V s {T n ) for any / G V s {T n ). In particular, Q e F G Vec&\ The operation 
Q e is a formalization for V s (T n ) of the usual finite-dimensional operation that allows one 
to consider the composition Xof of a vector X G l / ect < - s ' ) and diffeomorphism / G T> s (T n ) 
as a vector in Veci^. It denotes the shift of a vector, applied at the point f(x), to the 
point x with respect to global parallelism of the tangent bundle to torus. 

Lemma 3 (|5j) The following relations hold: 

TRg-i{Q g X) = Q e (TR g -iX); (21) 

and 

TRg(Q g -iX) = Q e (TR g X). (22) 
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Proof. By the above formulae we see that Q e (TR g -iX) sends the point m G T n to 
(m, X(g~ l (m))). On the other hand, Q g X = (g(m), X(m)), and hence 

TR r -i{Q g X) = (m,X(g-\m))) 

so that (|2ip is proved. Formula (1221) follows from (12T1) under the replacement of g by <? _1 . 
□ 

It should be pointed out that Q g is the global parallel translation at g on T> s (T n ) 
generated by global parallelism on T n . It turns out that Q g is the parallelism of Levi- 
Civita connection of metric (1T81) . Thus, for a smooth vector field Y(t) along a smooth 
curve g(t) in V s (T n ), the covariant derivative %Y(t) at a time instant t* is defined as 

jY(t) lt=t * = ^(Q^Yit))^. (23) 
As usual, a smooth curve g(t) in T> s (T n ) such that 

= 0, (24) 

is called geodesic. In the framework of Lagrangian approach to hydrodynamics PQ[2], 
describes the motion of so-called diffuse matter on T n without external force (the case 
with non-zero external force is constructed in analogy with that of perfect incompressible 
fluid below). For such g(t) introduce the vector v(t) G T e V s {T n ) (i.e., the .£F-vector field 
v(t, m) on T n ) by the formula v(t) = g{t) o g^it). It is shown that v(t) satisfies the Hopf 
equation (sometimes called Burgers equation without viscous term): 

d 

—v + {vV)v = 0. (25) 
dt 

Remark 2 It is shown in [2] that if g(t) is a geodesic, then for every f G T> s (T n ), the 
curve Rfgif) is also a geodesic. 

Lemma 4 Let g(t) satisfy ( 1241 and x G W 1 be an arbitrary specified vector. Then l x g{t) 
satisfies ( |24|) as well, where l x is introduced in Remark^ 

Proof. Note that the parallel translations given by operators Q g and Qf (see (120]) ). are 
commutative. By construction f t l x g(t) = f t (g(t,m) + x) = Qi x g{t) and Qi x (Qg(t*)9(t)) = 
Ql x g{t*)4jlx9{t)- Taking into account definition f l23l) of covariant derivative we obtain 

isU*) = Qu%m = o. □ 

Introduce the subspace (3f C TfV s (T n ) as TRf{3 with (3 introduced in ^2] (see, e.g., 
(fT3"|) ). Having done this at every / G V s (T n ), we obtain the smooth subbundle ft of 
TT> s (T n ) that in constructions below will be considered as constraint. This constraint is 
holonomic, i.e., the distribution ft is integrable. The integral manifold going through e 
is the submanifold and subgroup T> s ^(T n ) in T> s (T n ) that consists of i7 s -diffeomorphisms 
preserving the volume (see details in [2]). 

Notice that for Y G ftf the vector Q e Y may not belong to ft e = ft. 

Consider the map P : TT> s (T n ) — > ft determined for each / G T> s (T n ) by the formula 

P f = TR f o PoTRj 1 , 
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where P = P e : Vect {s) = T e V s {T n ) -> /? = (3 e = T e V s ^(T n ) is the orthogonal projection 
introduced in (1131) . It is obvious that P is D^(T n )-right-invariant. An important and 
rather complicated result (see [2]) states that P is C°°-smooth. 

It is a routine fact of differential geometry that the covariant derivative ^Y(t) of a 
vector field Y(t) along a curve g(t) in T> s p {F n ) is defined by the relation 

Let F(t,g,Y), Y £ T g V s ^{T n ), be a (force) vector field on D*(T n ). Consider a curve 
g(t) satisfying the equation 

^(i) = F(t,<7(t), <?(*))■ (26) 

Denote by u(t) the curve in T e T> s ^(T n ) (i.e., a divergence free vector field on T n ) obtained 
by right translations of vectors g(t), i.e., u(t) = g(t) o ^ _1 (t) = TR~^g{t). It is shown in 
[2] that u{t) satisfies the Euler equation 

ju + (u ■ V)u - gradp = TR~ l F{t, g(t),u(t, g(t))). (27) 

It should be pointed out that (1261) is Newton's second law with force F that describes 
the motion of perfect incompressible fluid on T n under the action of force 

TR^F(t,g(t),u(t,g(t))) 

depending of the "configuration of fluid" . Recall that a curve satisfying ( |26l) with F = 
is a geodesic. 

Remark 3 It is shown in [2] that if g(t) is a geodesic in T> s ^(T n ), then for every f 6 
T> s ^(T n ) the curve Rjg(t) is also a geodesic, but nothing like Lemma^ is valid in this 
case. 

If F is a right-invariant vector field on P^(T n ) such that F e = F, where F is a 
divergence free vector field on T n , then (J27J) turns into 

d 

—u + (u ■ S7)u — gradp = F. (28) 

Consider the map A : D s (T n )xf n TV s (T n ) such that A e : R" ->• (3 e = T e V s ^(T n ) C 
T e V s (T n ) is equal to A introduced by ([19]), and for every g E V s (T n ), the map A 9 : R™ -> 
T g T* s (T n ) is obtained from A e by means of the right translation, i.e., for X E R n : 

4(X) = Ti? 9 o A e (X) = (A o g)(X). (29) 

Since ^4 is C°°-smooth, it follows from w-lemma that A is C°°-smooth jointly in X E R n 
and g E V s (T n ). In particular, the restriction A : X>^(T n ) x R n -> TV^T 11 ) is C°°- 
smooth and the right-invariant vector field A(X) is C^-smooth on T> s ^(T n ) for every 
specified X E R n . 



8 



For any point m G T n , denote by exp m : T m T n — > T n the map that sends the vector 
X G T m T n into the point m+X in T n , where m+X is obtained modulo factorization with 
respect to the integral lattice, i.e., by the following procedure: We take a certain point in R n 
corresponding to m G T n , (denote it also by m) and X G R n = T m R ra , then we identify R n 
with T m R n = M. n , find m + X in 1" and pass from R n to T n by factorization with respect 
to Z n . The field of maps exp at all points generates the map exp : T e T> s (T n ) — > T> s (T n ) 
that sends the vector X G T e V s (T n ) (i.e., a vector field on T n ) to e + X G V s (T n ), where 
e + X is the diffeomorphism of T n of the form (e + X)(m) = m + X(m). 

Consider the composition exp o A e : R™ — > T> s {T n ). By the construction of ^4 e for any 
X G R" we get expo A e (X)(m) = m + X, i.e., the same vector X is added to every point 
m. Thus, obviously, expo A e (X) G T> s JT n ) and so expo A e sends R n to £>*(T n ). 

Let w(t) be a Wiener process in R n defined on a certain probability space (fi,jF, P). 
Introduce the process 

= exp o A e (aw(t)) (30) 

in P*(T n ), where cr > is a real constant. By construction, for u) G £1 the corresponding 

sample trajectory w!f\t) is the diffeomorphism of the form w!f \t)(m) = m + aw^t) 
so that the same sample trajectory aw^it) of is added to each point m G T n . In 
particular, this clarifies the fact that W^ a \t) takes values in T>*(T n ). 

Remark 4 Note that for a "specified u G fi" (^.e. ; a.s. /orw G f2j and for specified t G R 
£/ie vaZue of aw^it) is a constant vector in R n . T/ien, /or "given" uj and t, the action of 
W^\t) coincides with that of l a w u {t) ( see RemarkUl)- 

4 Viscous hydrodynamics 

In this section we take s > | + 2 so that the diffeomorphisms from V s (T n ) and T> s ^(T n ) 
are C 2 -smooth and Vect^ consists of C 2 -smooth vector fields. 

Everywhere below we use the same process constructed from a specified Wiener 
process w(t) in R n by formula (l30l . If, in the formula, there are several random elements 
with subscript uj, this means that they all are taken at "the same" wefl, i.e., sometimes 
the formula may be considered as description of a non-random element depending on the 
parameter uj G Q. 

Let g(t) be a solution of (124]) with initial conditions g(0) = e and g(0) = Vq G 
T e T> s (T n ). It is shown in [2] that such a solution exists on a certain time interval i G [0, T] 
(for the sake of convenience we take a closed interval inside the domain of g(t)). Recall 
that g(t) is a flow of diffuse matter without external forces. Consider v(t) = g(t)o o~ 1 (t) G 
T e T> s (T n ). This infinite-dimensional vector considered as a vector field on T n , will be also 
denoted by v(t,m). Recall that this vector field satisfies the Hopf equation ff25l) . 

Consider a process on V s (T n ) of the form rj(t) = W ia \t)og(t), t G [0,T], where W {cT) {t) 
is introduced in fl30|) . In finite-dimensional notation, 77 (t) is a random diffeomorphism 
of T n of the form rj(t,m) = g(t,m) + o~w(t) modulo the factorization with respect to 
integral lattice. Introduce the process £(t) = rj(T — t), i.e., in finite-dimensional notation, 
£(t, m) = g(T — t,m) + aw(T — t). 
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Since w(t) is a martingale with respect to its own "past", one can easily derive from 
the properties of conditional expectation that D*£(t) = g(T — t, m) = v(T — t, g(T — t, m)) 
and so PD*D£{t) = ^g(s)\ s=T _ t = 0. 

Consider the random process 

Us) = C(s) o r\t) = W&(T - s) o g (T - s) o g -\T - t) o (W&(T - t))~\ 

Notice that the random diffeomorphism (W^(t)) -1 acts by the rule 

{W^it))- 1 ^) = m-aw(t). 

Obviously £t(t) = e. The finite-dimensional description of this process can be given as 
follows. 

By construction 

m = C\t, m)) = g(T - t, C\t, m)) + aw(T - t). 

Then 

g(T - t, C\t, m))=m- aw(T - t) 

and so 

£-\t, m) = g-\T -t,m- aw(T - t)). 

Thus, 

&(s, m) = £(s, g' l {T -t,m- aw(T - t)) = 

g(T — s, g~ (T — t,m — aw(T — £))) + crw(T — s). 

Notice that £ t (i, m) = m — aw(T — t) + aw{T — t) = m, i.e., indeed, £ t (t) = e on V s (T n ). 

Since £ t (t) = e, the "present" a-algebra is trivial and so the conditional expecta- 
tion with respect to it coincides with ordinary mathematical expectation. Hence, using 
the relation between v(t) and g(t) and definition of D*, one can easily derive that 

DM*h=t = E(v(T -t,m- aw(T - t))) = E{Q e TR^ )(T _ t) v{T - t)) (31) 

(note that here t is specified and the derivative is taken with respect to s). 

Introduce on T n the vector field V(t,m) = E(y(t,m — aw(t))). We also denote this 
field, as an infinite dimensional vector, by V(t) = E(Q e TR^ a) ^v{fj). Formula (IBTj) 

means that D^ t {s)\ s =t = V{T — t). 

Theorem 1 The vector field V(T — t, m) satisfies the Burgers equation 

—V(T -t,m) + (V(T - t,m) ■ V)V{T -t,m)- — VV(T -t,m) = 0, (32) 

where V 2 is the Laplace- Beltrami operator which on the flat torus coincides with the 
ordinary Laplacian. 
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Proof. For t G [0,T] and u G Q, introduce the curve Ct,a>( s ) in s G [0,T] depending 
on parameter u, by the formula 

Ct,M = R^ ){T _ t) 9(T - s, g-\T - t)) = g(T - s, g~\T -t,m- aw(T))). 

Note the only difference between Ct,ui(s) and £t iU ,(s): there is the stochastic summand 
aw(T — s) in the expression for while it is absent in Q,uj(s). This means, in 

particular, that Ct,uj{s) is an a.s. smooth curve with random initial condition 

Note that £Ct,u,(s)\ s =t = -TR^ {a)( ^_ t v(T - s). Since g(T - s) is a geodesic, from 
Remark[2]it follows that "for almost all specified u" (i.e., a.s. for u) G Q) the curve Ct.^( s ) 
is also a geodesic, i.e., j~d(t,u{s) — 0. By Remark @] the action of the diffeomorphism 
w!f\t) coincides with that of l a w w (t)- Hence, by Lemma H] the curve 

is a.s. geodesic as well, i.e., ^£l(aw u (T-t))Ct,u( s ) = °- Note tnat 

Recall that EQ e TR~\ a) v(T -t) = V(T - t) and D4 t {s)\ s=t = V(T - £) (see above). 
Then from the above arguments and construction we derive that 

D*D£ t {s)\ s=t = D*V(T -t,Ct(s))\ s =t = ~E (J^^l(aw^(t))(tA s )\s=t^ = 0. 

But since D*£ t (s)\ s =t = V^(T — £), by formula (jHJ) the backward derivative D*V(T — t, 6( s )) 
coincides with the left-hand side of ( l32i) . Hence ( l32i) is fulfilled. □ 

Now, let us turn to the case of incompressible fluids. Let g(t) be a solution of (|26|) 
on X>*(T n ) with F = and with initial conditions g(0) = e and g(0) = u G T e £>*(T n ). 
As well as in the case of diffuse matter, it is shown in [2] that such a solution exists 
in a certain time interval £ G [0, T] (for the sake of convenience we again take a closed 
interval inside the domain of g(t)). Recall that g(t) is a flow of perfect incompressible 
fluid without external forces. Consider u(t) = g{t) o g^it) G T e V s ^{T n ). This infinite- 
dimensional vector considered as a divergence free vector field on T n , will be denoted 
u(t,m). Recall that this vector field satisfies the Euler equation f|28|) without external 
forces (see §3$). 

Since W {<J \t) takes values in V s ^(T n ) (see $3]), we can repeat on V^(T n ) the above 
constructions for V s (T n ), i.e., introduce rj(t) = W^'{t) o g(t), where t G [0,T], and 
= rj(T — t) (i.e., in finite-dimensional notation £(f, m) = g(T — t,m) + crw(T — 
t j). One can easily see that D*£(t) = g{T — t,m) = u(T — t,g(T — t, m)) and so 

PD*D4(t) = ig{s)\ 8= T-t = on X*(T*). 

As well as above the process £ t (s) = £(s) o £ _1 (t) has the property £ t (t) = e. Its 
finite-dimensional description is quite analogous to the case of V s (T n ). 

Introduce on T n the vector field U(t,m) = E(u(t,m — aw(t))) (a direct analog 
of V(t,m)). We also denote this field as an infinite dimensional vector by U(t) = 
E(Q e TR^ m u(t)). 
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Lemma 5 The vector field U (t, m) is divergence free. 

Proof. By construction, for an elementary event co G f2, (W^^)^) -1 is a shift of the 
entire torus by a constant vector. Hence, Q e applied to TR^ ((T) ^ u(t) means the parallel 
translation on torus of the entire divergence free vector field u(t) by the same constant 
vector back. Thus Q e TR^ (tT) ^u(t) is a random divergence free vector field on the torus. 
Hence its expectation is divergence free. □ 

So, U(t) G T e r>*(T n ). In particular, we have proved above that 

DM*)\s=t = U(T-t). (33) 

Since nothing like Lemma H] holds on V s ^{T n ), we have PD^D^ t (s)\ s=t = D*U{T — 
t^t{s))\ s =t 7^ (the value of this mean derivative is calculated in Remark |5] below). 
Hence there is no analogue of Theorem [TJ We can prove only the following 



Theorem 2 Vector field U(t,m) satisfies the following Reynolds type equation (see, e.g., 

m.j 



d_ 

m 



U + E 



(u ■ V)m ) (t, m — aw(t)) 



a 



V 2 U — gradp = 0. 



(34) 



Proof. It follows from Ito formula that 

du 



du(t, m — aw{t)) 



(t, m — o~w(t))dt + yV 2 u(t, m — aw{t))dt — an' dw(t), 



where V 2 , as well as above, is the Laplace-Beltrami operator and u! is the linear operator 
of derivative of u in m e T n . 

Recall that u(t,m) satisfies the Euler equation without external force, i.e., || = 
—P((u ■ V)u). Since 

d d d 

E(-u(t, m - aw(t))) = ^Eu(t, m - aw{t)) = —U(t) 



and E(o~(Vu)dw(t)) = 0, we derive that 







dt 



d 



—U = E[-u(t,m-aw(t)) 



dt 



E 



a 



-P[(u-V)u)(t , m — o~w(t)) + —\7 2 u(t, m — o~w(t)) 



-E 



(u ■ V)u J (t, m — aw(t)) 



V 2 U + gradp. 



So, (J33D is satisfied. □ 

There are usual methods for transforming ( l34l) into the standard Reynolds form (see 
|12j). For a divergence-free vector field X(m) on T n (i.e., for a vector X G T e T> s ^(T n )) 
introduce the random divergence free vector field 



U x (t, m) = X(m - aw(t)) - E(X(m - aw(t)) 
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(i.e., the vector 

U x (t) = Q e TR^ M(t) X - E{Q e TR^ {a){t) X) 
in T e V s ^{T n )). Evidently, for X = u(t), we obtain 

U u (t) (t, m) = u(t, m — aw(t)) — U (£, m) 

and so 

u(t, m — aw{t)) = U (t, m) + U u (t) (t, m) 
and EU u (t\(t,m) = 0. Then one can easily see that 

E([(u ■ V)u](t, m - aw(t))) = (U ■ V)U + E[(U u{t) ■ V)U u{t) }. 
Thus, ( 134)) transforms into 

q-U + {U- V)U - °— V 2 U - gradp = -E[{U u{t) ■ V)U u{t) } (35) 

which is the standard form of Reynolds equation. It differs from the Navier-Stokes type 
relation with viscosity ^- by the external force — E[(U u (f) • V) U u (t)\ that depends on u(t, m), 
not on U(t,m). Recall that fl35l) describes turbulent motion of fluid if the dependence of 
E[(U u (t) ' V)L r u ( 4 )] on U is given (say, derived from experimental data, see [12]). 

Remark 5 Note that for introduced above, formula ( |33l) reads that D^ t {s)\ s =t = 
U(T — t). Then, taking into account formula (JED, one can easily derive that 

PDM(s)s=t = PD*U(T - s, Us))s=t = ^U+{U- V)U - yV 2 f/ - gradp. 

Thus, eq. (I3"5"j) implies that PD 1f D 1f ^ t {s) s=t = —PE[(U u ^r-t) • ^)U u (T-t)]- 

Our next aim is to show that a slight modification of the above scheme of arguments 
allows us to annihilate the external force in (I3"5"|) by introducing a special random force 
field on V^(T n ) into (|26l). 

For a random divergence free a.s. H s+l - vector field X w (m) on T n (i.e., for a random 
vector X u G T e T>* +1 (T n ) C T e P*(T n )), construct the random vector field IJx ul {t,'m) 
which, for any uj G Q, is given by the formula 

U Xu (t,m) =Xu{m- o-Wu{t)) - E(X u) (m - aw^t)). 

Introduce the non-random H s vector field PE[(Ux u • V)C/x w ] an d then construct the 
random vector ^(t,X w ) in T e V s ^(T n ) by the formula 

= Q e TR w ^ (t) PE{(U Xuj ■ V)U Xu ). 

Note that PE[(Ux u • V)E/xJ and so $u(t, XS) lose the derivatives, i.e., they are iP-vector 
fields only since X u (and so UxS) is H s+1 - Thus dojitiX^) is well-posed only on an 
everywhere dense subset T e r , * +1 (T n ) in T e V s ^{T n ). 
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Now introduce the right-invariant force vector field g, Y u ), where Y u G T g V s ^{T n ), 
on T> s ^(T n ) that at g G X>* +1 (T n ) and G fi is determined by the formula 

&,(t, r w ) = TRg$ u {t, TRfYj, 

where TR~ X Y^ is a divergence free a.s. i/ s+1 -vector field. 
Consider the equation 

^9u,{t)=$ u (t,g u {t),g u {t)) (36) 
at 

on V s ^(T n ) whose right hand side is well-posed on the everywhere dense subset V^ +1 {T n ) 
in J) s (T n ). Note that (|36|) has no diffusion term and so it is an ordinary differential 
equation with parameter lo G Q. Here we do not investigate solvability of (136]) but 
suppose that for the initial condition g w (0) = e and <7 W (0) = u G T e "D^ +1 (T n ), it a.s. 
has a unique i7 s+1 -solution g u (t) which is a.s. well-posed on a non-random time interval 
t G [0,T] for a certain T > 0. Consider the divergence free a.s. i7 s+1 -vector field u w (t,m) 
on T n given by the relation g u (t) = u^it, g u {t)). The analog of above-mentioned vector 
U now takes the form 

U(t,m) = E(u w {t,m - <rw w {t))) = EQ e TRr\ a) u w (t). (37) 

As well as in Lemma [5] it is easy to see that vector field fl37j) is divergence free. 
Theorem 3 The divergence free vector field U jiuen by (1371) . satisfies the Navier-Stokes 

2 

equation without external force and with viscosity 

— U+ (U- V)U- ^-V 2 U-gradpi = 0. (38) 

Proof. Note that for the random field u' w (t, m) of linear operators and the random field 
u'^(t, m) of bilinear operators (the primes denote derivatives of u in m G T n ) the stochastic 
integrals f* u'^t, m)dw u (t) and f u'^(t,m)(dw u} (t),dw u; (t)) = J*tiu^dt = f* W 2 u u dt are 
well-posed. Then by applying standard arguments to the Taylor series expansion of u u , 
one can easily see that the Ito formula is well-posed for u u (t, m — aw^it)) and so 

E(du u (t, m - o-w u (t))) = E (^ u oj{t, m - aw u (t))dt + ^-V 2 ^^, m - aw u (t))dt \ . 

From (JSH]) it follows (see ([26]) and ([27])) that = -P[{u u ■ V)uJ[ + &,(*, Thus, 
in the same manner as in the proof of Theorem [2] and deriving ([35]) . we obtain 

— U(t, m) = E(—u w (t, m - aw u {t))) = 
ot at 



-E 



[v>u ■ V)0 (t, m - aw u (t)) 



2 

V 2 U + gradp 



a -72 



2 

+£Q e T/^\ CT)w ^(t,^(t)) = (39) 



2 

_(U • V)U + yV 2 U + gradp - E[{U Uui[t) ■ V)U uUt) ] 
+EQ e TR£ (cHt) $ u {t,u u (t)). 
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But by construction and by formulae (12T|) and (|22|) we get 

EQ e TR^ ){t) !S u {t,u u (t)) = 
EQ e TR^ ){t) Q e TR w ^ {t) PE[(U u ^ t) • V)^ (i) ] = 
EQeTR^^ TR^^Q^^PEliUu^t) ■ V)U u „ {t) ) = 

PE[(U u „ (t yV)U Uu>(t) }. (40) 

Recall that for any divergence free fields U, the vector fields J^-U and V 2 U are diver- 
gence free. Hence, gradjo in (139]) is taken from relation f|T4|) for E([(u u ■ V)u w ](t, m — aw^it))), 
i.e., 

PE{[{u^ ■ V)u u ]{t, m - awuif))) = E{[{u^ ■ V)itJ(t, m - aw w (t))) - gradp. 
Introduce gradpi and gradp2 by relations 

P(U • V)U = (U ■ V)U - gradpi 

and 

PE[(U uUt) ■ V)C/ Mw(t) ] = E[(U uUt) ■ V)U uUt) ] - gradp 2 . 

Evidently, gradp = gradpi + gradp2 (i-e., to within additive constants p = pi + P2)- Thus 
(138]) follows form ([39]) and (0OD in the natural form |U+ (U- V)U- 4v 2 U-gradpi = 0. 
□ 
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